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Abstract 

This work investigates the effect of finite-alphabet input constraint on the secrecy rate of a multi- 
antenna wiretap channel. Existing works have characterized maximum achievable secrecy rate or secrecy 
capacity for single and multiple antenna systems based on Gaussian source signals and secrecy code. 
Despite the impracticality of Gaussian sources, the compact closed-form expression of mutual infor- 
mation between linear channel Gaussian input and corresponding output has led to broad use of the 
Gaussian input assumption in physical secrecy analysis. For practical considerations, we study the effect 
of finite discrete-constellation on the achievable secrecy rate of multiple-antenna wire-tap channels. 
Our proposed precoding scheme converts the underlying multi-antenna system into a bank of parallel 
channels. Based on this precoding strategy, we develop a decentralized power allocation algorithm based 
on dual decomposition to maximize the achievable secrecy rate. In addition, we analyze the achievable 
secrecy rate for finite-alphabet inputs in low and high SNR regions. Our results demonstrate substantial 
difference in secrecy rate between systems given finite-alphabet inputs and systems with Gaussian inputs. 
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I. Introduction 

Wireless communications, with increasing coverage and applications, are vulnerable to po- 
tential security compromises such as passive eavesdropping and active jamming. Traditionally, 
network planners have relegated system security considerations to higher network layers of the 
OSI protocol stack through authentication and cryptography. However, in recent years, there have 
been growing research interests in the security analysis of wireless systems from a physical layer 
and information theoretic perspective. In a wiretap channel environment originally introduced 
by Wyner [?], a sender "Alice" wishes to transmit a secret message to the intended receiver 
"Bob" in the presence of a passive eavesdropper "Eve". Wyner [?] showed that when the Alice- 
to-Eve channel is degraded from the Alice-to-Bob channel, Alice can encode and send secure 
messages to the destination at a non-zero secrecy rate. In [?], a generalization for the non- 
degraded broadcast channel is proposed, and in [?], secrecy capacity of a Gaussian wiretap 
channel is shown to be achievable by adopting a random Gaussian codebook. In [?], the secrecy 
capacity of a multi-antenna Gaussian wiretap channel is shown to be achievable using a suitable 
input covariance matrix and by encoding the message using a Gaussian random codebook. 

For both single- and multi-antenna Gaussian wiretap channels, the codebook that achieves 
secrecy capacity turns out to be Gaussian. However, such codebooks are not implementable in 
practice. In real world systems, input codebook consists of finite set of equi-probable constellation 
points (e.g. M-QAM, M-PAM etc.). Therefore, in contrast to the Gaussian codebook, practical 
wiretap codes must consist of finite- alphabet symbols. Because of this constraint, the achievable 
secrecy rate for a finite- alphabet input scenario would differ from the secrecy rate achieved by 
a Gaussian codebook. 

A recent work [?] considered the effect of M-PAM input on the secrecy rate of a Gaussian 
wiretap channel and provided the necessary condition for power allocation to maximize the 
achievable secrecy rate. In [?], results were also extended to the case of parallel Gaussian wiretap 
channels. In [?], we have investigated the effect of finite- alphabet input on the ergodic secrecy 
rate of a multiple-input single-output and single-eavesdropper (MISOSE) system. To continue 
our progresses in this work, we investigate the effect of finite- alphabet input in a more general 
setting of a multiple-input multiple-output and multiple-eavesdropper (MIMOME) system. The 
specific contributions of the work are summarized below : 
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• In order to quantify the effect of finite- alphabet input on MIMOME systems, we propose the 
application of precoding matrix to transform the MIMOME channel into a bank of parallel 
channels. 

• We propose a power allocation optimization framework based on decentralized dual decom- 
position technique to maximize the achievable secrecy rate of MIMOME systems with an 
arbitrary but known input distribution. 

• We provide secrecy rate analysis of MIMOME systems with finite- alphabet inputs at low 
and high SNR regions. Our findings suggest that similar to the Gaussian wiretap channel, 
proper transmission power should be diverted at high signal-to-noise ratio (SNR) in case 
of finite- alphabet input albeit with different effect. 

We organize the rest of the paper as follows. We begin with the system model in section [IlJ In 
section [III] we propose a linear precoding scheme that transforms the MIMOME wiretap channel 



into a set of parallel channels. Based on this precoding scheme, we reformulate the secrecy 



rate problem for an arbitrary input distribution. In section IV we develop a decentralized power 
allocation algorithm based on dual decomposition that maximizes the achievable secrecy rate for 
an arbitrary distribution. In section [V] we further consider the special case of Gaussian input and 
present a modified water-filling power allocation strategy by considering the secrecy constraint. 
We then extend the modified water-filling power allocation scheme to analyze the secrecy rate 



for an arbitrary input distribution in both low and high SNR regions. In section VI we present 



numerical results before concluding with section VII 



II. Preliminaries and System Description 

Throughout this work, we use notations tr(.), det(.) and superscript {-} H , respectively, to 
denote the trace, the determinant, and the conjugate transpose of a matrix. 



A. System Model 

We consider a MIMO (multiple-input multiple-output) wiretap system model in which the 
transmitter (Alice), the intended receiver (Bob), and the passive eavesdropper (Eve), respectively, 
have m a (transmit), m b (receive), and m e (eavesdrop) antennas. Denote the received signals at 
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Bob and Eve as y b and y e , respectively. Their received signals are written as 

y b = H 6 x + n 6 (1) 
y e = H e x + n e (2) 

where H 6 G C mbXma and H e G C mjXroa denote, respectively, the flat-fading MIMO channels, 
from Alice-to-Bob and from Alice-to-Eve. The noise G C mb and n e G C me are zero-mean 
identity matrix variance complex Gaussian random vectors independent of each other. The data 
signal is x G C ma transmitted by Alice in the form of x = Ws, in which W is a linear 
precoding matrix. We denote s as a random vector with zero mean entries and identity correlation 
matrix. We constrain the total transmission power by a peak level Pt, i.e., tr{K x } < P T , where 
K x = E{xx H } is the covariance matrix of the transmitted signal vector. 

The secrecy capacity of the above system is achievable by using a Gaussian random codebook 
[?]. First, let A y B denote that A — B is non-negative definite. The secrecy capacity is the 
solution of the following optimization problem 

maximize logdet (I + H b K x Hf ) - logdet (I + H e K x Hf ) 

K x 

subject to : K x b 0, K x = (3) 
tr(K x ) <P T . 

In order to realize an achievable secrecy rate for an arbitrary input signaling, we will generalize 
the objective function. Instead of the optimal Gaussian signaling as used in the above optimization 
problem, we replace the objective function with the more general form of / (x; y&) — I (x; y e ), 
where / (x; y) represents the mutual information between the input and output vectors x and y. 

The optimization problem described in ([3]) is a non-convex optimization problem (except 
for the special case of m& = m e = 1 [?], [?]). Even for simple cases, the objective function 
possesses a number of local maxima. Therefore, the optimum value of K a . is in general not 
known. In particular, for an arbitrary input signaling without any closed-form mutual information, 
solving the above optimization problem can be a difficult task. Thus, instead of solving the 
above optimization problem, we consider a particular linear precoding scheme similar to the one 
provided in [?]. The proposed linear precoding scheme first transforms the MIMOME channel 
into a bank of parallel channels. This step allows us to gain a better understanding of the 
effect of finite- alphabet input on the secrecy of the system. It enables us to gain better insights 
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into future system implementation. Such linear precoding scheme is justifiable from practical 
system implementation perspective. In addition, at high SNR, such precoding scheme is known 
to achieve the capacity of a MIMOME system [?]. 



B. Preliminaries 

Definition 1: Similar to [?], we define the following subspaces 

S h = null(H b ) ± n null (H e ) 
S b , e = null(H 6 ) ± n null(H e ) ± 
S e = null(H 6 ) n null(H e ) ± 
S n = null(H 6 ) n null (H e ). 

□ 

In fact, subspace S b corresponds to the class of input with non-zero gain towards the direction 
of Bob only. Subspace S b ^ e corresponds to the class of input with non-zero gain in the direction 
of both Bob and Eve. S e corresponds to the class of input with non-zero gain in the direction 
of Eve only. Finally, <S n is the subspace with non-zero gain in the direction occupied by neither 
Bob nor Eve. Define k = rank ( Hf ^ and hence dim (S n ) = m a — k. In addition, 

we define, r = dim (Sb) and s = dim («S r , e ). Therefore, dim («S e ) = k — r — s. 

Definition 2: We recall the following definition of generalized singular value decomposition 
(GSVD) [?], [?] that we will use for our analysis. The GSVD of the pair (H 6 , H e ) takes the 
following form 



H ft = S h 



ma—k 



1 o ] *f 



H. = * P S f 



Q- 1 



H 



where * a G C maXma , * 6 G C m » xmb , and * e e C meXme are unitary matrices, fl G C kxk is a 
non-singular matrix. S fe G C mbXk and S e G C meXk have the following form 
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k— i — s 
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k—r—s 


I 








S e = s 





D e 





m e —k+r 












D 6 = diag ({&!, . . . , b s }) and D e = diag ({e a , . . . , e s }) are diagonal real matrices with dimension 
s x s. The diagonal entries of D 6 and D e are arranged in the following orders: 

o < h < . . . b s < 1 

1 > e 1 > . . . > e s > 

and 

b 2 i+e 2 i=l, for i = l,...,s. 

□ 

III. Secrecy Rate Problem Formulation 
A. Proposed Precoding Strategy 

Given the above two definitions in hand, we propose the following precoding matrix 



W = * a BP 1/2 (4) 

where, B is defined as follows 



k 


m,a-k 


ft 












P = diag ({pi, . . . ,Pm a }) is a diagonal power allocation matrix. The proposed precoding matrix 
in Eq. is similar to the precoding strategy defined in [?]. However, unlike the precoding 
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k—r- 



m a —k 



k—r- 



m e —k+r 





diag ({&i v /p fc _ r _ s+1 , . . . , b s y/Pk-r}) 

diag ({y/Pk-r+u y/Pk}) 

(10) 

k—r—s s m a —k+r 

diag ({y/pl, • ■ • , y/Pk-r-s}) 

diag ({e lv /^ 



(IT) 



scheme in [?], the diagonal elemens {p^ in the power allocaion matrix can have different 
values. As will be evident in the subsequent sections, such differentiation is extremely important 
due to the finite nautre of the input constellations considered in this work. 

By using the above precoding strategy, the system equations of ([!]), @ become 



y& 



*6 S b 



k m a — k 



I 



k m a —k 



I 



Kb 



P 1/2 s + n e . 



(6) 



(7) 



Pre-multiplying ([6]) and |7]) with and \l>^, respectively, we have the following equivalent 
equations 



H b s + n 6 
H e s + n e , 



(8) 
(9) 



where we use the notations y b = *f y&, y e = *f y e , fib = *f n b and n e = i&f n e . The new 



equivalent channel matrices H& = E& 



I 



I 



P 1 / 2 are specified in 



P 1 / 2 and H e = S e 
Eqs. (10j) and ([11}. 

From Eqs. ( fTO] ), ( [TT| ), we observe that the new system equations ([8]), (|9]) in fact transform 



the MIMOME system ([T]), (|2]) into a bank of parallel channels. Fig. [T] shows the resulting 
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parallel channel model. In this parallel channel model, input symbols si, . . . , Sk- r - s are only 
observed by Eve, symbols Sk- r -s+i, • • • > Sfc-r are received by both Bob and Eve, whereas symbols 
Sk-r+i, ■ ■ ■ ,Sk are only received by Bob. Finally, symbols Sk+i, ■ ■ ■ , s ma are lost by receivers of 
both Bob and Eve. 



B. Reformulation of Secrecy Rate Problem 

In this section, we relax the secrecy capacity problem of Eq. (|3]) using the precoding matrix 



presented in Sec. Ill- A We generalize the problem into an achievable secrecy rate problem for 
arbitrary input distribution. To this end, we present the following proposition. 

Proposition 1: Define I (j) = I (s; + n), v = k — r — s, p = {pi, . . . ,p ma }, = 
diag (fl H n) . When input s of a MIMOME system is a random vector with zero mean entries 
and identity correlation matrix, by using the precoding matrix W defined in Eq. Q, we can 
achieve the following secrecy rate for an arbitrary distribution of s 




subject to : ^ P»+i + ^2 Pj < p t 

i:bi>ei j=k—r+l 



Proof: See Appendix |A} ■ 
We observe that, the application of the proposed precoding matrix can transform the MIMOME 
problem into a distributed secrecy rate problem for a bank of parallel channels. We note that, 
for a given alphabet set, in general the above optimization problem should be jointly optimized 
over the input probability distribution and the power allocation. However, practical modulation 
constellations are generally constrained to be equi-probable. Therefore, here we will consider 
equi-probable input alphabet and focus on power allocation optimization. Next we propose a 
power allocation algorithm to solve the above optimization problem. 

IV. Power allocation algorithm for arbitrary input distribution 

Even though the above optimization problem is convex for Gaussian input, for an arbitrary 
input distribution, this is in general not the case. In addition, the lack of closed-form expression 
for mutual information makes the problem even more difficult to solve. In order to find an efficient 
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sub-optimal solution, we will revert to the decomposition technique. We note that, without the 



sum power constraint, the optimization problem ( fT2| ) can be decoupled into a number of parallel 
problems each involving only one variable pi. However, the sum power constraint compels us to 



solve a larger optimization problem jointly involving multiple variables. For problem such as ( |T2| ) 
involving a complex constraint, a dual decomposition method [?], [?] based on the Lagrangian 
of the objective function enables us to readily decompose the problem into a number of parallel 
sub-problems each involving a single variable. These subproblems are linked through a master 
problem that updates the dual variable during each iterations of the subproblems. By introducing 



a dual variable \l, we can write the Lagrangian of the optimization problem ( |12| ) by relaxing the 
coupling constraint as follows 



maximize 
p 



i.bi>e, 



1 



k 



Pv+i - 1 



el 



-l■ — r-L^ y J 



j=k—r+l 



Pj 



/' ! Y Pu+i + Y p j ) + t iPj 

\i:bi>ei j=k—r+l 



(13) 



subject to : Pi> 0. 

Note that, the above optimization problem is decoupled in terms of We obtain the following 
two subproblems. 

Subproblem 1 : for all i such that hi > e\ 

\ / p 2 M 

(14) 



maximize 

Pi>0 



I I '''' Pu+i ) I ( —— l>»+, 



u+i 



^Pu+i 



Subproblem 2 : for all j — k — r + 1, . . . , k 

maximize X 

Pi>0 



(15) 



Subproblems 1 and 2 are linked through the master problem which is the dual optimization 



problem of ([T2j). The master problem updates the value of the dual variables fi. Let p* denote 
the solution found from the subproblems. The master problem can be written as follows 
Master Problem : 



minimize >^ 

r.bi>ei 



b 2 



u+i 



X 



e 2 

% * 



P 



j=k—r+l 



(16) 



\i:bi>ei j=k-r+l 
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Subproblem 2 is a convex problem and can be solve optimally. However, in general subproblem 
1 is not convex. Even though subproblems 1 and 2 only involve single variable, in most situations, 
a closed form or analytic expression for the mutual information is not known for an arbitrary 
input distribution. In order to solve subproblems 1 and 2, we resort to a recent result on finite- 
alphabet research [?], [?] that relates the mutual information and the minimum mean square 
error (MMSE) at the receiver through 

^M = mmse(p). (17) 
dp 

Note that the function "mmse" for different discrete constellations (e.g., M-PSK, M-QAM etc., 



where M is the number of constellation points) has been given in [?]. Using Eq. ( [17] ), the 
optimum value of p, L in subproblems 2 can be solved from the following equations 
1 / 1 A 

mmse — p* J — p = 0, for j — k — r + 1, . . . , k. (18) 



Eq. ( [18] ) can be further expressed as follows 

p* = ojj mmse -1 (min {1, fiUj}) , for j = k — r + 1, . . . , k (19) 

where we used the fact that mmse _1 (l) = 0. 

For subproblem 1, we can derive the following necessary condition for optimality 

b 2 ( b 2 \ e 2 ( e 2 \ 

— mmse ^—pl +i — mmse ^—p* u+i - fi = 0, 

for i : bi > 6i (20) 



Next, we propose a sufficient condition for the optimality of the solution of Eq. ( [20] ). In this 

regard, we define the following MMSE difference function 

b 2 ( b 2 \ e 2 ( e 2 \ 

/mmseD(p, K ^u+i) = mm se P* u+l — mmse p* u+i . (21) 

\UJ u+ i J U>u+i \^>u+i J 

We note that, a similar condition has been proposed in [?] [see Theorem 6 in [?]] for a parallel 
Gaussian wiretap channel with M-PAM inputs. 

Proposition 2: If the MMSE difference function in Eq. ( [2T| ) admits a unique zero p' and is 
strictly monotonically decreasing for < p < p', then the optimal solution {pt +i } of Eq. ((20]) 
can be given as follows 

when b 2 < e 2 
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Pu+i 

when &? > e 2 



K 



p' if p = 
— — mmse I — — p* H 
Proof: See Appendix |Bj 



e 2 



mmse 



e? 



+ i 



fi, <p <p' if /x > 0. 



In following proposition, we describe a condition for the MMSE difference function to admit 
a unique zero solution. 

Proposition 3: If g(p) = pmmse(p) is a strictly uni-modal function and fej > ej, then the 
MMSE difference function f mmse D(p, e i; u F+! ) admits a unique zero solution for p > 0. 

Proof: See Appendix |Cj ■ 

In Figure [6j we present the plot of g(p) = pmmse(p) vs. p for BPSK, QPSK, 16-QAM, and 
64-QAM input constellations. We observe that g(p) of all four constellations shows strictly uni- 
modality in the region of interest. Furthermore, in Figures [7J [8} |9} and 10 we illustrate graphically 
the MMSE difference function f mmse B(p,a) vs. power p for various values of a = e}/b}. Here, 
we also assume uj v +i = 1 without any loss of generality. We observe that, the MMSE difference 
function exhibits strictly monotonically decreasing behavior in the range < p < p' for all 
four constellations. Based on Propositions [2] and [3j and Figures [6 10 we see that the optimality 
criterion described in proposition [2] holds for common constellations of BPSK, QPSK, 16-QAM, 



and 64-QAM. Therefore, the solution of subproblem 1 in Eq. ( [19] ) will be unique for at least 
these four constellations. 

The combined algorithm to solve the optimization problem ( [12] ) is presented in Algorithm ([T|). 



Algorithm 1 Dual decomposition algorithm for ( [12] ) 



1) Initialize dual variable p > 0. 



2) Find solution of subproblem 1 and 2 by solving ( [20] ) and ( [T9| ), respectively. 

3) Update dual variable as 

(k > 
i:bi>ei j=k-r+l 

4) Go to Step 2 until stopping criterion is reached. 



We note that the master problem in Eq. ( fT6[ ) is differentiable with respect to the dual variable 
p. Therefore, in Step 3 of the above algorithm, we use the gradient method to update the dual 
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variable. Here, the parameter a > denotes an appropriate stepsize, which can be either a 
constant or time-varying. In our simulation, we choose a fixed stepsize a. If the stepsize is 
sufficiently small, then the solution of the above algorithm will converge to the solution of the 
optimal dual variable /i*. A detail description on the choice of stepsize and stopping criterion 
can be found in [?], [?] and the references therein. If the original optimization problem in Eq. 
( [12] ) is convex, then the duality gap will be zero. Therefore, the solution of the dual problem in 
Eq. ([16]) will also provide the optimal solution of the original problem. However, if the original 



problem is not convex, then there exists a positive duality gap and the solution of Eq. ( [16] ) will 
be a suboptimal solution of the original problem. 



V. Secrecy Rate Analysis of MIMOME System 
A. Gaussian Input Case 

In this section, we will present the power allocation problem for the special case of Gaussian 
input distribution based on the framework considered above. Even though the result for Gaussian 
input is well known [?], [?], [?], results in this section will provide additional insight for the 
finite- alphabet input scenarios to be considered later. 

For Gaussian input, the MMSE equation is simply 

mmse(7) = (22) 



1 + 7 



Let us denote {p\} as the optimum power allocation for Gaussian input. Based on Eq. ( [22] ), 
solutions of subproblem 1 can be found by solving the following equation 



b 2 e 2 



1 



(pI + i) 2 + 



b l ~ e ; Uv+i 
Thus, the optimal power allocation is 

0. 



b 2 _ e jPu+l 



0. 



(23) 



p. 



v+i 




+ 4 



Wi 



u v +i 



ifi< 
ifi> 



u, 



v+i 



el 



V K - e\ 



Similarly, based on Eq. (22), solution of subproblem 2 can be found as 



Pi 







0. 


if — < Ui 




n 


1 




— Wi, 


if — > Ui 







(24) 



(25) 
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We notice that, the solution of subproblem 2 presented in Eq. ( f25| ) can also be obtained from 
Eq. (|2~3]) for by replacing b 2 with 1 and ef with 0. In the absence of eavesdropper, the solution of 



this simplified problem is the famous water-filling solution as given in Eq. ( |25| ), where the base 

level is and water level is the inverse of the dual variable — . When the security constraint 

fi 

is present, however, the problem takes on an interesting structure. In such case, we can still 



1 



consider the water level as — . However, we will use 

fi 



V+l 



as a base level to take into account 
the additional constraint due to secrecy. For the parallel channels j = k — r + 1, . . . , k, with 



b 2 - e? 



components only towards Bob's direction, the base level will still be Uj (since, b 



f_„ = 1 and 



-J-U 



0). Similar to water-filling, we will allocate power only when the water level is above 



the base level. However, the power level in this case will not be the difference between the water 
level and base level, i.e. 



b 2 - e? 



Instead, it will be a non-linear function of the difference 



as given in Eq. p4| ). As a result, the achievable secrecy rate can be written as 

\ 



m 



i:bi>ei 



1 1 + — p 



0J. 



u+i 



1 + 



0J. 



■ p 



v+i 



j=k— r+1 



P, 



(26) 



B. Low SNR Approximation 

1) Second order optimal signaling: For second order optimal signaling [?], the first and the 
second order derivatives of the mutual information achieved at zero SNR matches with those 
achieved using Gaussian input. In general, quadratic symmetric signaling such as QPSK or any 
other signaling distribution that can be written as a mixture of QPSK (i.e., M-QAM, for M > 4) 
are second order optimal. For second order optimal signaling, the low SNR approximation of 
MMSE (i.e., the first derivative of mutual information) is the same as that of Gaussian signaling. 
Hence, we can use the same water-filling power allocation solution as presented in Eq. ([24]), 



in Section V-A 



2) Non-second order optimal signaling: 1-D signaling schemes such as BPSK and M-PAM 
are not second order optimal. A low SNR approximation of such signaling is given in [?] as 



mmse (p) = 1 - 2 p + o (p 2 ) . (27) 

Based on the above equation, low SNR power allocation {p' ow } for non-second order optimal 
signaling are given below in two cases: 
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1 ^V-t-i 

0, when - < 



"Tr - A* To 2 when - > 



0, when — < Ui 

P X r={ Ut / 1 \ ? , i = *-r + l,...,*. (29) 

— /i Wi , when - > Ui 

2 \p J fi 

The achievable secrecy rate at low SNR can be approximated as follows 

R l r= J2 (b*-e*)(p%-( P ^) 2 ) + £ (pr-M ow ) 2 )- (30) 

i:bi>ei j=k—r+l 

We notice that, similar to the case of Gaussian and second-order optimal signaling, we can 

use a water-filling strategy with water level — and base level { ^ vJr ' 1 anc [ assuming 6f = 1 

H bf-ef J 

and e 2 -_ v = for j — k — r + 1, . . . , k. The power level is still function of 

1 w v+ i 



» b\ - ej 

which is the difference between water level and base level. The functions are now slightly 



different in (Eq. (|28) and (|29j) 



C. High SNR Approximation 



1) rank(H e ) = m a : In this case dim(S&) = r = 0. For the power constraint in Eq. ( fT2] ), 
we get the following complementary slackness condition 



^2 ^u+iPu+i -Pt) =0. (31) 



\i:bi>ei 



Based on proposition [2j the optimal power p* for maximum secrecy satisfies p* < p'. Therefore, 
at very high SNR when Pt — > oo, secrecy rate for finite- alphabet input is maximized by using a 
fraction of the total available power. The power constraint inequality becomes a strict inequality. 
In the above complementary slackness condition, we attain p = 0. Denoting high SNR power 



allocation as pf gh , we re-write Eq. f20|) as 



b't mmse ( b\ pJJ ) = e\ mmse ( e\ pjg ) . (32) 
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In [?], based on the sub-optimum estimator s(y,p) = axgmin \y — y/josj the following 
MMSE approximation at high SNR is found: 



2 



f d 

mmse(p) « K exp < — ^ p i» , (33) 

in which K is a constant and d is the minimum distance between two signaling points in the 
discrete unit variance input constellation. 

Reference [?] also provided a table containing formula for calculating d of different finite- 



alphabet constellations. Using Eq. ( |33] >, we obtain the following high SNR approximation of 
power allocation 



Pv+i = d 2 ^tf^ e 2\ > fOT i SUCh that b i> e i- (34) 



4 uJu+i 

In [?], it was observed that, at high SNR regime, the power allocation for parallel Gaussian 
channel with finite-alphabet demonstrates a channel inversion characteristic. In other words, 
stronger channels receive less power allocation. This is because the mutual information of a M- 
ary constellation cannot exceed log 2 M bits/s/Hz, Thus, there is little incentive to allocate more 
power to a channel once the mutual information is near saturation. Instead, additional power is 



better allocated to weaker channels for higher rate. In Eq. (34), we observe a similar channel 

(b 2 - e 2 ) 

inversion phenomenon, although in this case, the effective channel — — . This is in sharp 

contrast to both the water-filling power allocation at low SNR regime and the power allocation 



for Gaussian input (see Eq. ([28])), where the power allocation was proportional to the effective 
channel. 

2) rank (H. e ) < m a : Based on Eq. ( [33] ), mmse(p) decays exponentially to zero as p — > oo. 
From Eq. ([18]), therefore, we find that p — > as P T — > oo. Hence, for the subset of parallel 



channels, % : bi > e^, Eq. (34) will still provide a high SNR approximation of power allocation. 



The subset of parallel channels j = k — r+l, . . . ,k are in the subspace <S&. For these channels, 
there are no components in Eve's subspace and a similar channel inversion style power allocation 
can be achieved as presented in [?] with an effective channel — . For an M-ary constellation, 
mutual information for these channels will become close to log 2 M at high SNR. 

In summary, for the subset of parallel channels i : bi > a channel inversion based power 

b 2 - e 2 

allocation based on the effective channel — will be performed, whereas for the set of 

parallel channels j = k — r + 1, . . . , k, a channel inversion type power allocation based on 
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effective channel gains l/ooj will be performed. A high SNR approximation of the achievable 
secrecy rate for the case rank(H e ) < m a is given by 

fl«*«rlogM + \ Z (— PIS) 



1,2 \ / 2 

hioh \ „ / Cj high 



i:6i>ei 

VI. Numerical Results 



Pv 



(35) 



A. Transmitter with Accurate Eavesdropper CSI 

Without any loss of generality, we assume equal noise power level at receivers of both Bob 
and Eve. We also assume that Alice has full CSI of both Bob and Eve. Our numerical results 
average over 500 channel realization, where each entry of both Bob's and Eve's channel matrices 
is i.i.d. complex random Gaussian variable with zero mean and unit variance. 

In Fig. [2[ we present numerical test results for a 5 x 5 x 5 MIMOME system. In this test, we 
ensure that each realization of H e is non-singular, i.e., rank(H e ) = 5 = m a . Since, rank(H e ) = 
m a , there is no parallel channel only directed at Bob. In addition to result for the proposed power 
allocation (PA) algorithm, we also present results obtained using the water-filling PA of Section 



V-A and results from equal power over all channels (uniform PA) [?]. We also present high SNR 
approximation results for every tested constellation as well as the low SNR approximation result 
for BPSK. For other constellations (QPSK, 16-QAM and 64-QAM), low SNR approximation 
gives the same result as the water-filling PA. 

As seen in Fig. [2j at high SNR, power allocations according to water-filling and uniform 
strategies would drop the secrecy rate to almost zero. This result is intuitive. For finite-alphabet, 
the achievable mutual information at high power approaches the saturation value of log 2 M. 
For Gaussian input, however, the mutual information or the capacity increases monotonically 
with increasing power. Since both water-filling and uniform strategies assume a Gaussian input 
distribution, at high SNR both schemes would use more power to transmit signals. This strategy 
drops the secrecy rate asymptotically to zero as the difference in mutual information between 
Alice-to-Bob and Alice-to-Eve narrows with at very high SNR. Fig. [2] also indicates that the high 
SNR approximation analysis closely matches the secrecy rate at high SNR regime. Similarly, the 
secrecy rate at low SNR is also closely approximated by the low-SNR approximation analytical 
result. 

Fig. [3] presents test results for a 5 x 5 x 3 MIMOME system. Because in this case rank(H e ) < 
m a , there are parallel channel components only directed towards Bob. Therefore, results obtained 
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in this case are different from Fig. |2j Specifically, we observe that, even though secrecy rates 
using water- filling and uniform PA schemes at high SNR drops from the maximum attainable 
value, they do not approaches zero as in Fig. [2} The reason is due to the existence of channels 
directed only at Bob. For the parallel channels with components along both at Bob's and Eve's 
direction, the difference in mutual information between Alice-to-Bob and Alice-to-Eve narrows 
at high SNR. However, for each of the parallel channels with components only in Bob's direction, 
the secrecy rate would approach log 2 M. Therefore, at high SNR the total secrecy rate approaches 
r log2M. Moreover, our proposed PA algorithm can achieve some additional non-zero secrecy 
rate from the channels directed at both Bob and Eve, in addition to those parallel channels only 
in Bob's direction. 

Both Fig. [2] and [3] indicate that when using Gaussian assumption for finite-alphabet inputs, 
there is a threhold SNR above which the achievable secrecy rate starting to decrease. The value of 
this threshold SNR is higher for the higher order constellation. Therefore, by adaptively switching 
to the next higher order modulation format once we reach the threshold point for a particular 
constellation, we will be able to use Gaussian water-filling PA algorithm even for finite- alphabet 
inputs. We also notice similar threshold SNR for uniform PA. In addition, uniform PA also 
achieves secrecy rate close to the water-filling PA. Therefore, a very simple but near optimal 
strategy would be to use uniform power allocation and start switching to the next higher order 
modulation once above each threshold SNR. 

B. Transmitter with Partial Eavesdropper CSI 

Thus far, our analysis assumes that Alice possesses full channel information of Eve. In practice, 
however, Eve's CSI or even the presence of a passive Eve is difficult to determine. Therefore, 
in this section, we will numerically evaluate the scenario when Alice only has access to partial 
(statistical) information regarding Eve's channel state. In particular, let Eve's channel consists 
of 

H e = H e + E e . 

Here, H e is Eve's mean CSI known to Alice whereas E e is the CSI uncertainty which is modeled 
as zero mean white Gaussian noise with variance <j\, i.e., E e ~ CM (0, of I). 

In this case, if Alice performs power allocation based on the known observation H e by 
disregarding the uncertainty, the power allocation may not be optimal. In fact, Alice may even 
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lose secrecy. In the following proposition, we present the achievable ergodic secrecy rate for a 
given power allocation p. 

Proposition 4: For a given power allocation p based on the known observation H e without 
considering the uncertainty, the achievable ergodic secrecy rate -R sec (p) can be given as follows 

6? 



i? S ec(p) = Yl 

i:bi>ei . 



'Pv+i 



Here, e, ~ CM (0, a\ cUi), i = 1, . . . , k are i.i.d. random. Also, a\ 
Proof: See Appendix |D| 




(36) 



From Eq. ( |36| ) we notice that, data symbols Sfc_ r+ i, . . . , Sk - previously only observed by Bob, 
are also seen by Eve now. In addition, for the symbols s v+ i, . . . , s u+s - observed by both Bob 
and Eve, an additional uncertainty term e v+i is added to the mutual information of Eve. 

In Fig. [4]and[5j we present test results of achievable ergodic secrecy rate for different values 
of variance of channel uncertainty, for 5 x 5 x 5 and 5x5x3 MISOME systems, respectively. 
Here, Alice is only aware of the mean CSI H e and use this CSI for power control without 
considering the uncertainty component E e . We present test results for different values of available 
transmission power. Both Fig. [4] and [5] indicate that the achievable ergodic secrecy rate decreases 
with larger channel uncertainty. In addition, we observe from Fig. [4] that the achievable ergodic 
secrecy rate plots for transmission power P T > 5 stays the same. This is because Alice does not 
require all available power to transmit at higher SNR. Hence, the power allocation stays the same 
at higher SNR. However, in Fig. |5J we observe that the achievable ergodic secrecy rate in fact 
decreases with larger transmission power for P T > 10 dB. In this case, rank (H e j < m a . When 
Alice only uses this information for power control without considering the uncertainty, Alice will 
likely allocate more power to the bank of parallel channels which only has components toward 
Bob. At high SNR, Alice will allocate more powers to these channels. However, as shown in 



Eq. ( |36| ), because of channel uncertainty E e , Eve now also possesses components along these 
channels. In other words, Eve can also receive signals from these channels. As SNR grows large, 
the mutual information for finite alphabet saturates to log 2 M. As a result, Eve can receive nearly 
full data information in these channel and consequently, the achievable ergodic secrecy rate will 
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decrease at high SNR. 

We note that, in order to maximize the achievable ergodic secrecy rate under channel uncer- 
tainty, one needs to solve the following optimization problem 



maximize 
p 




i \ _r_/ q 



+ £ 

j=k— r+1 
s k 

subject to: J^2Vh + PJ - P t- 

i=l j=k—r+l 

Investigation of this optimization problem is beyond the scope of this work. However, in [?], a 
power allocation algorithm to maximize the achievable ergodic secrecy rate under arbitrary input 
distribution is provided for a multiple-input single-output and single-eavesdropper scenario. An 
extension of the power allocation algorithm provided in [?] can also be used to solve the above 
optimization problem. 

VII. Conclusion 

This work considers the effect of practical finite-alphabet inputs on the secrecy performance of 
an MIMOME system. Our investigation led to the application of a precoding matrix to convert 
the MIMOME system into a bank of parallel channels so as to reformulate the achievable 
secrecy rate problem. We proposed a decentralized dual decomposition and a corresponding 
power allocation algorithm to maximize the achievable secrecy rate based on the proposed 
precoding for channel transformation. We analyzed the Gaussian input as a special case and 
provided a water-filling inspired power allocation strategy. Furthermore, we derived analytical 
results of achievable secrecy rate based on approximations at low and high SNR scenarios. Our 
results show that power allocation strategy based on Gaussian input is far from optimal when 
applied blindly in finite- alphabet input situation and may in fact be very risky by driving the 
secrecy rate to zero at higher SNR. 

Appendix A 
Proof of Proposition Q] 

An achievable secrecy rate for an arbitrary distribution of input x can be obtained by solving 
the following optimization problem 
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maximize I (x; y 6 ) - / (x; y e ) 

K x 

subject to : K x y 0, K x = ( 37 ) 
tr(K x )<P T 

When we apply the precoding matrix W from Eq. the input of the system becomes x = Ws. 
The objective function of the above optimization problem can be written as follows 



I (x; y&) - I (x; y e ) 

-I (s; H b Ws + n 6 ) - / (s; H e Ws + n e ) 
=/ ( s; H fe s + h b ) - I ( s; H e s + n 



k—r—s 



Z(pi) + {^i Pk-r-s+i) 



( 6 i Pk-r-s+i) + Yl X ( p i 
i=l j=k—r+l 

k s 
j=k— r+l 



i=l 
k—i — s 



(38) 
(39) 
(40) 

(41) 
(42) 



Note that (40) follows from (39) since linear unitary transformation of channel outputs preserves 



mutual information. Because (J40j) represents the difference in mutual information between a 



subset of parallel channels, we can rewrite ( |40[ ) into the summation form of ( |4Tj ). 

Since s is a random vector with identity correlation matrix, K x = E [xx H ] = WW H , we 



have tr(K x ) = J2i=i u iPi- Therefore, we can reformulate the optimization problem in ((37]) as 
follows 



maximize 

{Pi} 



k—r—s 



i Pk—r—s+i 



) ~ I (e 2 i p k -r-s+i)] - Y ^(P*) + ^(P*) 



i=l 



i=k— r+l 



(43) 



subject to : 



i=l 



Note that the above optimization problem is a distributed power allocation problem for the 
bank of parallel channels shown in Fig. [TJ where the variables {jOj}'s are coupled through 
the total power constraint. Because Z(.) > 0, we have optimum power allocation p* = for 
i — 1, . . . , k — r — s. In addition, since 1(7) are monotonically increasing in 7, X (6| pfc_ r _ s+i ) — 
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X{efpk-r-s+i) — whenever &j < hence p* = 0, for all i with 5; < ej. Furthermore, from 
Fig. [TJ we observe that symbols Sk+i, • • • , s ma are transmitted towards the direction of S n and is 



not part of the optimization problem (43). Hence, it will be wasteful to expend any transmission 



power to transmit these symbols. Therefore, p* = for i = k+1, . . . , m a . By change of variables, 



p' = uiiPi, and then replacing p' with p, we obtain the optimization problem in Eq. (12) 



Appendix B 
Proof of Proposition [2] 

If bf < ef, the difference in mutual information between Alice-to-Bob and Alice-to-Eve 
channel will be negative. Hence, no secrecy is possible. As a result, the optimal power allocation 
should be zero. 

When bf > ef, it is possible to achieve a positive secrecy rate. If, in the solution of the 



original problem ( |12| ), the sum power constraint inequality becomes a strict inequality, then 



p = according to the KKT condition [?]. Such condition can occur at high SNR regime. 



Therefore, the solution of problem (20) will be achieved by the power allocation p' that renders 



the MMSE difference function in Eq. pi) zero. If the MMSE difference function admits a 



unique zero solution, then the solution of Eq. (20) will be optimal, i.e., p* u+i = p'. 



If on the other hand, the sum power constraint in the original problem admits an equality, 
then fj, > 0. If the MMSE difference function is monotonically decreasing for < p < p', then 
the function / mmse D(j>, e i; oo v+i ) — /i will be zero for a power allocation value p* that is less 
than p' . Now, if MMSE difference function is strictly monotonically decreasing, then the power 



allocation solution will be unique. Therefore, the solution of Eq. ( |20| ) will be optimal 



Appendix C 
Proof of Proposition [3] 

Define, p\ = b- L p and p 2 = eip. First, we will prove that the difference function, gu(pi, P2) = 
g(pi) — g(p2) admits a unique zero for p > 0. 

Assume that the strictly unimodal function g(p) is strictly monotonically increasing for p < m 
and strictly monotonically decreasing for p > m. When, p 2 < p\ < m, the difference function 
firj(pi) P2) cannot be zero due to the strictly monotonically increasing nature of the function 
g(p). Similarly, when p 1 > p 2 > m, the difference function gu(pi, P2) cannot be zero due to the 
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strictly monotonically decreasing nature of the function g(p). Therefore, the difference function 
9d(pi,P2) can be zero only when p 2 < m < p\. 

Suppose, for p = p', the difference function gu(hp, eip') equals zero. For p > p', g(bip) > 
gibip') and g[tip) < g(e,ip') due to the strictly unimodal properties of g(p). Hence, when p 2 < 
m < pi, the difference function <7d(pi,P2) = g{pi) — g{p2) cannot be zero for p > p'. The same 
argument holds for any point p < p' . Therefore, when 6, > ej, the difference function gn(pi, P2) 
admits a unique zero for p > 0. Hence the MMSE difference function / mmse D(p ; e», u u+i ) has 
a unique zero solution for p > 



Appendix D 
Proof of Proposition [4] 

With H e = H e + E e , Eq. @ can be written as 

y fi = H R s 



*?E P Ws 



n f 



(44) 



Here, H f 



*fH e W 



I 



P 1 / 2 . Now, the second term in the above equation can be 
expressed as *f E e Ws = E e PV 2 s , where, E e = E e B and E e = *f E e * a . We can decompose 
E e as a sum of two matrices, where the first matrix is a diagonal matrix containing the diagonal 
elements of E e . The second matrix contains the non-diagonal entries of E e and contains all zero 
elements in its diagonal. The decomposition is given below 



e x 

e 2 







.. 


.. 


•. 


"• 


• • h 


.. 


.. 


.. 


.. 


.. 











+ 





e 2 i 



ei2 




e/ci efc2 
e(fc+i)i e(fc+i) 2 



ei fc ... 

g 2fc ••. 












e me fc ... 



(45) 



Pi 

Here, ~ CM (0, cr^j). Based on the above decomposition and replacing with — , the 



i-th entry of the vector y e can be expressed as follows 

1 



PiSi + gi + n e 



(46) 
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Here, g t = J~] 



en 



pisi and e'j can be expressed as 



1, i = 1, . . . ,k — r — s 

ej_„, i = v + 1, . . . ,u + s and v = k — r — s (4V) 
0, i = k — r + 1, . . . , m e 
Assume that the eavesdropper is performing conventional decoding on each of these parallel 
branch by considering the cross terms as a part of noise. Then the achievable ergodic secrecy rate 
for a given power allocation p can be written as in Eq. ( |36] ). If in case, Eve employs advanced 
decoding scheme, e.g. successive interference cancellation etc. then Eq. ( [36] ) will serve as a 
upper bound on the achievable ergodic secrecy rate. 
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Fig. 1. Precoding matrix W = ^aBP 1 ' 2 converts MIMOME channel to a bank of parallel channels 
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Fig. 2. 5x5x5 MIMOME system 
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Fig. 3. 5x5x3 MIMOME system 
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Fig. 4. Secrecy Rate of a 5 x 5 x 5 MIMOME system with partial Eve's CSI 




Fig. 5. Secrecy Rate of a 5 x 5 x 3 MIMOME system with partial Eve's CSI 
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Fig. 6. p vs, pmmseO) plot for BPSK, QPSK, 16-QAM and 64-QAM constellations. 
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(a) 




(b) 




(c) 



Fig. 7. MMSE difference function /d{p, a) for different values of a — ef/bf for BPSK 
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Fig. 8. MMSE difference function /d(p, a) for different values of a — ef/bf for QPSK 
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Fig. 9. MMSE difference function fjj (p, a) for different values of a — 
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Fig. 10. MMSE difference function fo(p, a) for different values of a = ef/b? for 64-QAM 
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